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Recently, Matsukawa et al. have discovered a new superconductor Pr2Ba4Cu70i5_<s in which 
metallic CuO double chains are responsible for the superconductivity. To investigate the su- 
perconductivity, we employ the d-p double chain model where the tight-binding parameters are 
determined so as to fit the LDA band structure. On the basis of the Tomonaga-Luttinger liquid 
theory, we obtain the phase diagram including the superconducting phase in the weak coupling 
limit. We also calculate the Luttinger liquid parameter K p as a function of the electron density 
n by using the Hartree-Fock approximation. With increasing n from quarter filling, K p in- 
creases, and then exceeds 1/2 when the superconducting correlation becomes most dominant. 
K p has a maximum at an optimal density between quarter- and half-filling. These results are 
consistent with the experimental observation. 
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Low-dimensional strongly correlated electron systems 
have attracted much interest due to the possible rele- 
vance to the high-X'c superconductivity. Although the 
two-dimensional Cu02 planes play essential roles for the 
superconductivity, one-dimensional (ID) CuO chains in- 
cluded in some families of high-X'c cuprates have also pro- 
vided an interesting testing ground for anomalous metal- 
lic states due to effects of the electron-correlation. 1-3 

PrBa2Cu40s (Prl24) is an excellent system to study 
the properties of ID CuO chains, because the electronic 
conduction in Cu02 plane is suppressed due to the so- 
called Fehrenbacher-Rice state formed by the strong hy- 
bridization between Pr 4/ and O 2p orbitals. 4 By study- 
ing the anisotropy in the resistivity of a single crystal, 
it was clarified that the metallic conductivity is caused 
by conduction in CuO double chains. 5 Some researchers 
suggested that the CuO double chain might be in the 
Tomonaga-Luttinger liquid state. 6 

Pr 2 Ba4Cu70i5_,5 (Pr247) has both CuO double chains 
and CuO single chains and shows metallic conductivity 
at low temperatures owing to the metallic conduction in 
the CuO double chains. Furthermore, the carrier density 
can be varied by controlling the amount of oxygen defi- 
ciency in the CuO single chains. In the recent reports, 
the resistivity in Pr247 was investigated by changing the 
oxygen content, and superconductivity with T c ~ 15K 
was found after a reduction treatment. 7 ' 8 The NQR ex- 
periment revealed that the superconductivity is realized 
at the CuO double chains. 9 These experiments suggest 
that the material seems to show the possibility of novel 
ID superconductivity. 

Many theoretical works have been performed on ID 
strongly correlated electron systems confined to a dou- 
ble chain such as the two-chain model and the Ladder 
model. 10-21 At half-filling, the system can be described 
by a Heisenberg model whose ground state is a spin liq- 
uid insulator with a finite gap in spin excitation. 22 Away 
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from half-filling, the system becomes a metallic state 
which maintains a spin gap. 10-12 This behavior is ex- 
plained by the existence of electron pairs produced by the 
dominant fluctuations of the 4fcp charge density wave or 
the intcrchain-paring fluctuations. The paring fluctua- 
tion is expected to dominate over the other fluctuations 
in the weak coupling regime of the Hubbard model as 
well as in the strong coupling regime of the t — J model. 

Most of the theoretical works on the double chain 
systems have been done by using simplified single-band 
models such as the Hubbard model or the t — J model. 
As for the single chain systems, two-band models such 
as the d-p model, simulating a CuO linear chain, have 
also been extensively investigated. 23-28 They show that 
the charge fluctuation between d and p orbitals, which 
is not taken into account in the single-band models, has 
crucial effect on the superconductivity. In the case of 
CuO double chain (see Fig. 1), in addition to such ef- 
fect, the number of the Fermi points and the values of 
the corresponding Fermi velocity play crucial roles for 
the superconductivity. They directly depend on param- 
eters included in the d-p double chain model. Therefore, 
we need theoretical studies on the basis of a realistic d-p 
double chain model. 

In the present paper, we investigate the d-p double 




Fig. 1. Schematic diagram of Cu(3d x 2 _ y 2) orbitals and 0(2p CT ) 
orbitals in a CuO double chain. 
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Fig. 2. A typical energy dispersion relation for the CuO double 
chain. Solid lines are the tight-binding result (k) with t pd = 
1.6eV, t pp = 0.43eV, t dd = 0.12eV and A = 4.1eV. Closed 
circles are the LDA result for Y124. 34 Dotted line indicates the 
Fermi level with n = 0.57. 



chain model, simulating a CuO double chain composed 
of Cu(3d x 2_ y 2) orbitals and 0(2p a ) orbitals as shown in 
Fig. 1. The tight-binding parameters, i.e. the hopping 
integrals and the charge transfer energy are estimated to 
fit the LDA band structure. By using the weak coupling- 
theory, 29-31 we obtain the phase diagram including the 
superconducting phase. 

We also calculate the Luttinger liquid parameter 
K p 32 ' 33 as a function of the electron density for a finite 
Coulomb interaction between d electrons U d by within 
the Hartree-Fock (HF) approximation. 

The Hamiltonian of the d-p double chain model is 
given by 

H = t pd ^2(pl a d l+la + h.c.) + U d ^2 ndi-\n dll 

i.a i 

+ tpp J2(pLPi+i<y + h - c -) + e pYl 

i.a i,a 

+ t dd ^2(d f i(J d i+1<J + h.c.) + e d ^2d\ a d la , (1) 

where d\ a and p\ stand for creation operators for a elec- 
tron with spin a in a Cu(3d x 2_ y 2) orbital at site i and 
for a hole with spin a in a 0(2p a ) orbital at site i, re- 
spectively, and h d i a — d\ a di a . Here, t pd is the hopping 
energy between the nearest-neighbor d and p sites and t pp 
(tdd) is the hopping energy between the nearest-neighbor 
p (d) sites. The atomic energy levels for p and d orbitals 
are given by e p and e d , respectively. The charge-transfer 
energy A is defined as A = e d — e p . U d is the on-site 
Coulomb interaction between d electrons. 

In the noninteracting case with U d — 0, the Hamilto- 
nian eq. (1) is easily diagonalized. It yields a disper- 
sion relation: E ± (k) — \{t d + e p + 2(t dd + t pp ) cos A: ± 
^(A + 2{t dd - t pp ) cosfc) 2 + 16{t pd cos (fc/2)) 2 }, where 
k is a wave vector and E + (k) (E~(k)) represents the 
upper (lower) band energy. A typical energy dispersion 
relation £ ,± (fc) is depicted in Fig. 2. As the charge trans- 



fer energy is positive, A = e d — e p > 0, the upper band 
E + (k) corresponds to the d band, while the lower band 
E~(k) corresponds to the p band. 

We estimate the tight-binding parameters in the d-p 
model eq. (1) so as to fit the energy dispersion E ± (k) 
to the LDA band structure. Since the LDA results for 
Pr247 have not been available so far, we employ the LDA 
result for YaBa 2 Cu 4 8 (Y124) 34 in which CuO double 
chains are included with the same lattice structure as 
those in Pr247. 7 Therefore, the energy dispersions cor- 
responding to the CuO double chain for the both com- 
pounds are expected to be almost equivalent. Using 
the least squares method, we obtain the parameters as 
t pd ~ 1.6eV, t pp ~ 0.43eV, t dd ~ 0.12eV and A ~ 4.1eV, 
respectively. Here, the fitting weight of the d band is 
taken to be one hundred times as large as that of the p 
band. This choice gives better fitting near E-p than the 
fitting with equal weight for both bands. 

These estimated values are of the same orders of the 
well known values for the corresponding energies in the 
Cu02 plane. 35 The hopping energies satisfy the rela- 
tions: t pd > tpp > t dd > 0, which are consistent with the 
values of the overlap integrals between atomic orbitals 
(see Fig. 1). Then, we expect that the obtained tight- 
binding parameters are realistic to describe the electronic 
state of the CuO double chain in Pr247. 

Since the d band dispersion relation E + (k) has a dou- 
ble minimum (maximum) structure as shown in Fig. 2, 
the number of the Fermi points depends on the electron 
density n in the d band of the CuO double chain. There 
are two Fermi points for small n, while, four Fermi points 
for large n. This difference plays crucial role for the su- 
perconductivity as will be discussed later. 

The LDA calculation for Y124 34 has predicted that 
there are four Fermi points in the d band of the CuO 
double chain with the electron density n ~ 0.57. As for 
the case with Pr247, it is considered to be nearly quarter 
filling n ss 0.5 due to the Fehrenbachcr-Rice effect men- 
tioned before. 4 Then, the CuO double chain of Pr247 
is in the boundary between the two Fermi point system 
and the marginal four Fermi point system as seen in Fig. 
2 where the Fermi level for n — 0.57 is plotted. 

On the basis of the Tomonaga-Luttinger liquid the- 
ory, 29-33 various types of correlation functions show 
power-low dependence with critical exponents. These 
exponents are determined by a single parameter K p 
in the model which is isotropic in spin space. For 
the single-band model with two Fermi points, ±fc_F, 
the SC correlation function decays as ~ r k p , 
while the CDW and SDW correlation functions decay 
as ~ cos(2kpr)r^ i ' 1+K f'K When the system is in the 
Tomonaga-Luttinger regime, both of the charge and spin 
excitations are gapless (we label this regime as clsl). In 
this case, the SC correlation is dominant for K p > 1, 
while the CDW or SDW correlation is dominant for 
K p < 1. 

In contrast to the single-band model, the situation 
is rather complicated in the two-band model with four 
Fermi points, ifcpi and ikp 2 . When the ratio of the two 
Fermi velocities \vf 1 /vf 2 \ lS smaller than a critical value 




Fig. 3. Phase diagram on the n-t pp plane in the weak coupling 
limit for t pd = 1.6eV, A = 4.1eV and t dd /t pp = 0.28. Closed 
circle stands for the parameter point corresponding to the LDA 
calculation for Y124. 



Fig. 4. The Luttinger parameter K p as a function of the electron 
density n for several U d for t pd = 1.6eV, t pp = 0.43eV, t dd = 
0.12eV and A = 4.1eV. Inset shows K p as a function of U d for 
n = 0.57. 



~ 8.6, the low-energy excitations are given by a single 
gapless charge mode with a gapped spin mode (labeled 
as cls0). 29 ~ 31 In this case, the SC and the CDW correla- 
tions decay as ~ r 2K p and ~ cos[2(kp 2 — kp 1 )r}r~ 2K '' , 
respectively, while the SDW correlation decays expo- 
nentially. Hence, the SC correlation is dominant for 
K p > 0.5, while, the CDW correlation is dominant for 
K p < 0.5. When the ratio \vf 1 /vf 2 \ is larger than the 
critical value ~ 8.6, the low energy excitations are given 
by two gapless charge modes and two gapless spin modes 
(labeled as c2s2). 

In the non-interacting case with Ud = 0, the Luttinger 
liquid parameter K p is always unity. In the weak cou- 
pling limit Ud — > 0, i.e., in the limit K p — > 1, the SC cor- 
relation is most dominant in the clsO region (SC phase), 
while, the CDW or SDW correlation is most dominant 
in the clsl and c2s2 regions (normal phases). 29-31 In 
Fig. 3, we plot the phase diagram in the weak coupling 
limit on the n-t pp plane. The values of t p d = 1.6eV and 
A = 4.1eV together with the ratio tdd/tpp = 0.28 are 
fixed to the corresponding values in Fig. 2. 

As mentioned before, the CuO double chain of Pr247 
is expected to be in the boundary region between the 
clsl and the clsO phases. When the system is initially 
in the clsl phase, the electron doping effect brings about 
the phase transition from the clsl to the clsO (see Fig. 
3). This is consistent with the experimental observation 
in Pr247, where the superconductivity is caused by the 
oxygen reduction corresponding to the electron doping 
effect in the CuO double chain. ' 8 

When the distance between the two chains of a CuO 
double chain decreases, the hopping terms t pp and tdd 
are expected to increase. Therefore, the pressure effect 
might lead to the phase transition from the clsO to the 
clsl as shown in Fig. 3. This is again consistent with 
the recent experiment in Pr247 under high pressure, 36 
where the superconductivity is suppressed and finally 
disappears due to the pressure effect. 



Now, we calculate the Luttinger liquid parameter 
K p for finite Coulomb interaction Ud- Based on the 
Tomonaga-Luttinger liquid theory, K p is related to the 
charge susceptibility \c and the Drude weight D by the 
following equation: 33 



Kp = \^ X M 12 



with 



d 2 E a 



(2) 



(3) 



dn 2 ' " " d<t> 2 
where E g is the ground state energy per unit cell as a 
function of the electron density n and the magnetic flux 
4>. To calculate K p in the d-p model eq.(l), we estimate 
Eg within the HF approximation and substitute it into 
eq.(3). The obtained value of K p is valid up to the first 
order of Ud and is expected to be a good approximation 
in the weak and intermediate coupling regimes as shown 
previous works. 13 ' 14 ' 27 ' 28 

Fig. 4 shows K p as a function of n for several values of 
Ud/t p d with t p d = 1.6eV, t pp = 0.43eV, tdd = 0.12eV and 
A| n= o.57 = 4.1eV, where A is the renormalized charge- 
transfer energy within the HF approximation and is ex- 
plicitly given by A = A + As A depends on n and 
Ud, we determine the bare charge transfer energy A so as 
to fix A to the corresponding value in Fig. 2: A = 4.1eV 
at n = 0.57 for each Ud- 

As shown in Fig. 4, the clsl phase appears for n ^ 0.4, 
while the clsO phase appears for n 0.4. In clsO phase, 
Kp increases with increasing n and then have a maxi- 
mum at an optimal electron density n ~ 0.85. In the 
clsO phase with Ud/t p d ~ 4, the value of K p exceeds 
1/2 for a region of n around 0.85, where the SC correla- 
tion becomes most dominant as compared with the other 
correlations (SC phase). On the other hand, in the clsl 
phase, Kp is always smaller than unity and the CDW or 
SDW correlation is most dominant (normal phase). As 
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shown in the inset in Fig. 4, K p decreases monotoni- 
cally with increasing Ud\ the similar behavior has been 
observed in the {/-dependence of K p in the Hubbard lad- 
der model. 13 ' 14 

Experimentally, the superconductivity of Pr247 is ob- 
served only for a range of the oxygen reduction rate 
5 ^ 0.3 and T c shows a maximum at an optimal reduction 
rate S ~ 0.45. 8 When we assume the electron density of 
the CuO double chain in Pr247 with S = to be n ~ 0.4, 
we can guess n ~ 0.4+5 for finite 5. Then, for Ud/t p d = 4 
corresponding to Ud = 6.4eV, the SC region (n > 0.7) 
with K p > 1/2 as well as the optimal electron density 
(n ~ 0.85) with a maximum in K p seems to be consistent 
with the experimental observation mentioned above. 

To summarize, we have investigated the superconduc- 
tivity in the d-p double chain model, simulating a CuO 
double chain of Pr247, where the tight-binding parame- 
ters are determined so as to fit the LDA band structure. 
On the basis of the Tomonaga-Luttinger liquid theory, 
we have obtained the phase diagram in the weak cou- 
pling limit, which provides us a clear understanding of 
the doping and the pressure dependence of the supercon- 
ductivity in Pr247. 8 ' 36 

The Luttinger liquid parameter K p has also been ob- 
tained as a function of the electron density n using the 
Luttinger-liquid relation combined with the HF approx- 
imation. The doping dependence of K p is in good agree- 
ment with that of T c in Pr247, 8 when we assume that 
T c is finite for K p > 1/2 in the clsO phase and is mono- 
tonically increasing function of K p . Although the finite 
value of T c is not obtained in purely one-dimensional 
systems, an approximate value of T c could be estimated 
as a function of K p by taking into account of a three- 
dimensionality due to the effect of couplings between the 
double chains. We will report it in a subsequent paper. 
Even in the normal state above T c , a kind of transition 
between the clsl phase and the clsO phase is expected. 
In fact, such a transition has been observed in the doping 
and/or the pressure dependence of transport properties 
such as the resistivity and the Hall coefficient. 8 ' 36 In ad- 
dition to the transport properties, a spin gap is expected 
to exist in the clsO phase. A spin gap like behavior has 
been observed in the recent NQR experiment, 9,37 where 
(TiT) -1 is suppressed in the superconducting sample 
as compared with the non-superconducting sample even 
above T c . 

Taking account of the LDA calculation and our re- 
sult, we can expect that the CuO double chain of Y124 
is also in the clsO phase and shows the superconduc- 
tivity as well as Pr247. We note that a rapid de- 
crease in the resistivity below 15K has been observed in 
YBa 2 (Cui_ :E Zn a ;)408 with high Zn doping 0.05 < x < 
0.1. 38 In this case, the superconductivity in the Cu02 
plane is considered to disappear, but the CuO double 
chain might preserve the superconductivity. To be more 



conclusive, we need further investigation in this regard. 
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